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ABSTRACT 

The  experimental  results  for  convective  heat  transfer  and 
flow  friction  characteristics  of  some  screen,  plate-fin  and 
skewed  fin  matrices  of  various  geometry  are  presented,  The  heat 
transfer  data  were  obtained  using  the  transient  technique. 

The  results,  which  are  presented  in  graphic  and  tabular 
form,  should  be  of  practical  use  for  compact  heat  exchanger 
design. 

The  results  of  a  digital  computer  solution  giving  the 
maximum  slope  of  the  generalized  heating  curve  as  a  unique 
function  of  NTU  are  presented  in  tabular  form  in  Appendix  II, 

The  author  wishes  to  thank  Professor  C,  P.  Howard  of  the 
U,  S.  Naval  Postgraduate  School  for  his  encouragement  and 
guidance,  and  Lt,  J,  A,  Inglis,  RCN,  for  solving  the  Bessel's 
function  and  theoretical  generalized  heating  curve  equation. 
The  assistance  of  J,  S,  Beck  and  K,  Mothersell  in  fabricating 
the  equipment  used  is  also  acknowledged. 
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INTRODUCTION 

There  are  two  courses  open  to  the  designer  of  a  compact 
heat  exchanger  for  the  transfer  of  a  fixed  heat  load  consistent 
with  an  allowable  pressure  drop  utilizing  a  new  type  of  matrix 
arrangement  that  is  not  covered  in  present  heat  transfer  litera- 
ture. A  matrix  will  be  defined  as  a  porous  solid  possessing  a 
high  ratio  of  void  surface  area  to  bulk  volume,  and  a  void  geo- 
metry which  permits  the  flow  of  gas  or  a  liquid  through  it. 

The  first  course,  a  trial  and  error  technique,  consists  of 
analytically  predicting  the  heat  transfer  and  friction  character- 
istics of  the  matrix,  designing  and  testing  a  scale  version  of  the 
heat  exchanger,  and  modifying  the  design  until  the  required  perfor- 
mance is  achieved.  The  time  and  cost  involved  in  utilizing  this 
method  is  generally  unacceptable  to  the  designer. 

The  second  method  consists  of  generating  the  required  heat 
transfer  and  friction  data  on  a  small  module,  and  utilizing  the 
results  for  the  full  scale  design.  This  method  is  quite  flexible, 
since  many  different  matrix  arrangements  can  be  tested  at  mini- 
mum cost,  and  an  optimum  surface  chosen. 

One  of  the  satisfactory  methods  of  measuring  the  convective 
heat  transfer  characteristics  of  porous  solids  utilizes  the  so- 
called  transient  test  technique.  This  consists  of  introducing  a 
step  function  in  the  temperature  of  the  fluid  passing  through  the 
matrix,  which  is  initially  at  uniform  temperature,  and  recording 
the  matrix  temperature-time  history.  Using  the  results  of  the 
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analysis  developed  by  Gr.  L.  Locke  (l)*  for  the  determination 
of  NTU  from  the  generalized  heating  curve  of  maximum  tempera- 
ture gradient  versus  time,  it  is  then  possible  to  compute  the 
Stanton  number.  The  'Fanning'  friction  factor  is  obtained  from 
the  pressure  drop  across  the  matrix. 

The  objectives  of  this  thesis  were: 

a.  to  design,  fabricate  and  evaluate  the  performance 
of  a  test  section  to  utilize  the  transient  method  for  deter- 
mining heat  transfer  coefficients  of  small  modules,  and 

b,  to  obtain  the  heat  transfer  coefficients  and  fluid 
frictional  characteristics  of  screen,  triangular  fin  and  plate- 
fin  matrices  of  various  geometry. 


^Numbers  in  parentheses  refer  to  Bibliography  on  page  2! 
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EXPERIMENTAL  APPARATUS 

A  schematic  diagram  of  the  apparatus  is  shown  in  Fig8  8. 
'Basically,  the  system  consists  of  a  plastic  section  with  a 
sliding  drawer  for  holding  the  test  matrix,  a  flow  measuring 
system,  a  pressure  measuring  system,  a  compressor  to  provide  flow 
through  the  test  section,  an  air  heater,  a  temperature  measuring 
system  and  piping. 
Flow  Hate  Control 

The  compressor  is  run  at  a  constant  rate  of  approxi- 
mately 300  c.f.m.  A  sliding-gate  valve  on  a  T-section  of  the 
inlet  piping  of  the  compressor  is  used  for  coarse  control  of 
air  flow  through  the  test  section.  There  is  also  a  fine  control 
valve  for  making  minor  adjustments  to  the  air  flow  rate.  A  thin 
plate  orifice  section  for  flow  measurement,  designed  according  to 
A.S.M.E.  specifications  (3)»  is  included  in  the  inlet  piping 0 
Matrix  Heating  System 

a  small  blower,  with  a  heating  element  controlled  by  a 
115-volt,  60-cycle  AC  '^ariac',  is  positioned  to  heat  the  matrix 
when  the  drawer  is  removed  from  the  plastic  test  sections  A 
transitional  section,  with  flow  stabilizing  screens,  is  attached 
to  the  blower  to  accomplish  uniform  heating  of  the  matrix3 
Temperature  Pleasuring  System 

a  five -thermocouple  upstream- temperature  probe  is  located 
in  the  center  of  the  pipe  approximately  midway  between  the  ori- 
fice and  test  section.  The  upstream-downstream  air  temperature 


difference  is  measured  by  these  thermocouples,  connected  in 
series  so  as  to  oppose  the  5  center  thermocouples  of  a  13- 
thermocouple  tray  located  immediately  after  the  matrix. 

The  thirteen-thermocouple  tray  may  be  used  indepen- 
dently to  find  the  temperature  distribution  of  the  downstream 
air.  All  thermocouples  used  are  30-gauge  iron-constantan* 
The  temperature  difference  and  distribution  is  recorded  on  a 
Leeds  £c  Northrup  'Speedomax5,  Model  G,  emf  recorder  with  a 
variable  range  of  from  1-20  mv. 

Ambient  air  temperature  is  read  on  a  mercury  thermo- 
meter located  at  the  inlet  of  the  entrance  piping. 
Sliding  Drawer  Matrix  Holder 

The  sliding  drawer  matrix  holder  is  constructed  of 
half-inch  thick  plastic .  It  is  lined  with  balsa  wood,  al= 
lowing  for  a  3»08-inch  diameter  passage  for  holding  the  screen 
matrix,  when  testing  the  triangular  fin  matrices,  a  3.125- 
inch  square  passage  is  provided  for  holding  the  matrix. 
Pressure  Measurement  System 

The  pressure  drop  across  the  orifice  is  measured  on  a 
lJj-inch  i-tube  marometer  filled  with  water.  The  pressure  drop 
across  the  matrix  is  measured  on  either  a  50-inch  differential 
manometer  or  an  inclined  160  mm.  manometer  Cilled  with  water. 
Piping 

Ambient  air  is  drawn  through  a  3."'d-inch  diamet er  dpe 
to  the  test  section.  When  testing  the  triangular  fin  matrices, 


a  transition  piece  to  a  3. 125-inch  square  was  inserted  before  and 
after  the  test  sections  Galvanized  ducting  is  used  between  the 
test  section  outlet  and  compressor  inlet. 


METHOD 

A.  Summary  of  Theory 

The  theoretic-il  analysis  of  the  transient  heat  transfer 
behavior  of  a  porous  medium  was  first  presented  by  Shumann  in 
1929.  5humann  considered  a  homogeneous  porous  medium  at  a  uni- 
form temperature  through  which  a  fluid  of  the  sane  temperature 
is  flowing.  At  a  certain  instant,  the  temperature  of  the  enter" 
ing  fluid  is  assumed  to  change  to  a  higher  or  lorer  value.  The 
problem  is  then  to  find  the  temperature  of  the  fluid  and  solid 
as  functions  of  time  and  position  in  the  matrix. 

The  following  idealizations  and  initial  and  boundary 
conditions  have  to  be  satisfied  in  order  to  simplify  the 
problem  t 
Idealizations 

1.  The  fluid  specific  heat  and  viscosity  are  constant . 

2.  The  thermal  conductivity  of  both  fluid  and  solid  is 
zero  in  the  direction  of  flow  and  is  infinite  within 
the  fluid  and  solid  in  the  direction  normal  to  the  flow, 

3.  The  flow  is  steady. 

h»     The  porous  solid  is  honogenous9 
Initial  and  Boundary  Conditions 

1.  At   time  Q  -  0,  the  temperature  of  the  entering  fluid 
changes  instantaneously. 

2.  Initially,  the  core  is  at  a  uniform  temperature. 

3.  Mo  heat  passes  the  core  boundaries. 


Shumann  found  the  solutions  for  the  fluid  and  metal 


temperatures  to  be: 
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Shurnann' s  solution  requires  that  a  theoretical  curve  be 
generated  which  will  match  the  experimental  curve  to  deter- 
mine the  heat  transfer  coefficient, 

In  19!?0,  Locke  (l)  developed  a  method  in  which  only 
the  maximum  slope  of  the  theoretical  and  experimental  curves 
need  comparison.  This  method  eliminates  the  requirement  of 
a  large  number  of  theoretical  curves  and  the  possibility  of 
error  due  to  the  displacement  of  the  experimental  heating 
curve  either  vertically  or  horizontally,  Locke  derived  an 
expression  for  the  slope  of  the  generalized  heating  curve 
for  the  fluid,  (equation  1  above),  finding  that  at  x=L  for 
which  Z  becomes  equal  to  the  NTU: 


-(Nnitr; 


From  this,  it  was  found  that  the  maximum  slope  of  the  generalized 
heating  curve  is  a  unique  function  of  MTU.  Thus,  only  the  maximum. 
slope  of  the  experimental  heating  curve  need  be  found  in  order  to 
obtain  the  NTU  and,  thereby,  the  Stanton  number 8  This  method  has 
been  used  by  several  investigators,  and  their  results  have  been 

compiled  by  Kays  &  London  (2). 

i  it- 
Locke,  in  his  solution  of  Equation  3>  resorted  to  an 

approximation  of  the  modified  Bessel  function  which  in  the  low 

range  of  the  argument  was  subject  to  large  error.  A  solution 

to  Equation  3  has  been  carried  out  by  Lt.  J.  A.  inglis,  RON, 

using  the  exact  infinite  series  of  the  modified  Bessel  function, 

on  a  CDC  160U  digital  computer.  The  results  of  this  solution 

are  given  in  Appendix  II  and  were  the  ones  used  by  the  author 

in  reducing  his  experimental  data. 

B.  Experimental  Techhique 

Ambient  air  is  drawn  through  the  inlet  pipe  and  allowed 
to  flow  through  the  test  matrix,  while  the  flow  rate  is  ad» 
justed.  After  constant  flow  is  established,  the  sliding  drawer 
containing  the  matrix  is  pulled  out  of  the  test  section  and  into 
the  heated  air  stream. 

When  the  matrix  temperature  is  steady,  as  indicated  by 
monitoring  th*.  thermocouples,  the  drawer  is  rapidly  inserted 
into  the  ambient  air  stream  and  the  expei  iment  is  started., 

The  ambient  air  temperature,  orifice  pressure  drop  and 
matrix  pressure  drop  are  recorded,  and  the  transient  upstream 
to  downstream  temperature  difference  is  recorded  on  the  Leeds 
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St   Northrup  ' Speedomax1  Recorder.  After  the  temperature  differ- 
ence  has  returned  to  zero,  the  run  is  terminated. 
The  folio  win,,  data  are  recorded: 

a.  atmospheric  pressure 

b.  atmospheric  temperature 

c.  upstream  to  downstream  temperature  difference 

d.  pressure  drop  across  the  orifice 

e.  pressure  drop  across  the  matrix 
fo  orifice  diameter 

In  this  investigation,  the  temperature  difference  between 
the  inlet  fluid  and  the  downstream  fluid  temperature  was  recorded 
directly  so  the  maximum  slope  could  be  determined  from  the  ex- 
perimental trace  with  no  reolotting.  For  a  complete  set  of 
sample  calculations,  see  Appendix  I, 

C.  Presentation  of  iiesults 

To  avoid  the  confusion  often  encountered  when  a  large 
number  of  arbitrarily  defined  parameters  is  used,  the  same 
definitions  will  be  applied  throughout  to  the  various  matrix 
geometries. 

The  basic  data  and  results  for  each  surface  are  pre- 
sented in  both  tabular  (Tables  I  -  VI)  and  graphic (Figs.  l-£) 
forms,  as: 

NStKPr2/3  =  *2<!***) 
f  -  j)9(%e) 


The  Reynolds  number  is  based  on  a  hydraulic  diameter, 

defined  as  follows : 

Ur  G 
NKe=  IJL 

/< 
with 

^rh  =  Dh  -  UAcVA,  ft., 

G  =  wi'Ac9  IV  m  hr-ft  2 

The  Stanton  number  is  defined  as: 

Nq4  -  JL.  -  NTU  kJk 

Jl   Gcf       G 

where 

NTU  =   /wfcjp  -  number  of  heat  transfer  units  of  the 
matrix . 

The  friction  factor  is  defined  on  the  basis  of  an  equivalent 
shear  force  in  the  flow  direction  per  unit  of  heat  transfer  area. 
Entrance  and  exit  losses  are  included  in  the  friction  factor, 

f  »  AP  matrix  2gc  Ac/G-  VfA 

Actually,  Prandtl  number  was  not  a  test  variable,  as  the 
working  fluid  was  in  all  cases  air  at  moderate  temperatures, 
Qve"  a  moderate  range  of  Prandtl  numbers,  the  effect  of  Prandtl 
number  to  the  2/3  power  is  used  as  an  approximation, 

D.   Limitations  and  uncertainties 

Care  must  be  used  to  approach  the  idealizations  and  initial 
and  boundary  conditions  of  the  theoretical  analysis,  which  are 
restated  below: 
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Idealizations 

lo  The  fluid  specific  heat  and  viscosity  are  constant. 

2o  The  thermal  conductivity  of  both  fluid  and  solid  is 
zero  in  the  direction  of  flow  and  is  infinite  within 
the  fluid  and  solid  in  the  direction  normal  to  the  flow. 

3s  The  flow  is  steady. 

h»     The  porous  solid  is  homogenous. 

Since  the  fluid  used  is  air,  the  first  idealization  can  be 
approached  if  temperature  differences  are  small,  ks   the  maximum 
temperature  difference  is  approximately  20°  F.,  the  arithmetic 
average  will  be  plus  or  minus  1.0$  in  viscosity  and  negligible  in 
specific  heat. 

As  for  the  second  idealisation,  it  is  recognized  that  ther- 
mal conduction  along  the  flow  length  occurs  and  increases  as  the 
flow  rate  is  lowered.   Its  effect  can  be  lowered  by  shortening, 
or  dividing  the  length  of  matrix  in  the  flow  direction,  Itondt  (?) 
indicates  that  as  the  Reynolds  number  is  decreased  the  conduction 
effect  will  be  small  below  an  NTU  of  11,79,  have  no  effect  at 
NT J  ■  11,79,  and  increase  appreciably  as  the  Reynolds  number  is 
further  reduced. 

In  the  direction  normal  to  the  flow,  the  conduction  is 
negligible  sincg  the  temperature  gradient  at  any  section  is 
essentially  zero.  This  gives  a  condition  which  is  equivalent  to 
infinite  conductivity. 

The  third  idealization  is  approached  by  running  the  com- 
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pressor  at  rated  capacity  for  all  runs,  so  that  pulsing  does  not 
occur. 

The  fourth  idealization  is  approached  by  careful  selection 
of  the  porous  solid  to  be  tested. 
initial  and  Boundary  Conditions 

1.  at  time  9=0,  the  temperature  of  the  fluid  flowing 
through  the  matrix  changes  instantaneously. 

2.  Initially,  the  core  is  at  a  uniform  temperature. 
3 9  No  heat  passes  the  core  boundaries. 

The  first  condition  can  be  approached  by  rapidly  intro~ 
ducing  the  heated  matrix' into  the  flow  stream  of  constant- 
temperature  ambient  air» 

The  second  condition  is  satisfied  by  placing  screens  in 

>ter  duct  to  produce  a  uniform  velocity  profile  of  the 
heated  air.  The  downstream  side  of  the  matrix  is  monitored 
by  a  thirteen- thermocouple  tray  to  check  for  a  uniform  tem- 
perature distribution  ov.      section. 

The  third  condition  can  only  be  minimized  by  careful 
design  of  the  matrix  holder.  To  lower  the  transfer  of  heat 
from  the  matrix  holder  to  tl  e  matrix  on  cooling,  the  matrix 
holder  is  lined  with  balsa  wood. 

The  uncertainties  of  physical  constants  in  the  matrix 
and  duct  measurements  are  listed  below: 
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Quantity  Uncertainty  a1    20 si  odds 

cs  ±1.0* 

NPr  £1,2# 

/M  -1.0* 

^s  negligible 

L  &.# 

Ac  -1.0* 

A  -1.0* 
Instrumentation  uncertainties  follow: 

quantity  Uncertainty  at  20 ?1  odds 

ti-tfl  to .$% 

Aporifice  ~lo0* 

Apmatrix  ^2.0* 
Following  is  a  summary  of  the  uncertainty  analysis  of 
the  Reynolds  number,  friction  factor,  MTU  and  Stanton  number: 
i'or  Kepeatability 

Quantity  Uncertainty  at  20 si  odds 

%e  *2.0* 

f  -3.;^ 

NTU  i8e0-10' 

NSt  fe.O^lQ 

%tNPr2/3  iS.O-lO; 
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The  NTU,  No+  and  N-^Npj,2^  are  represented  with  a 
variable  range  of  uncertainties  because  NTU  does  not  vsry 
linearly  with  the  maximum  slope  of  the  cooling  curve, 


NTU 

Uncei 

•tainty  at  20:1  odds 

6 

Jic? 

10 

i  9.0% 

16 

m   3.056 

20 

-  9,0% 

for  Reproducability 

Quantity 

Uncertainty  at  20  fl  odds 

%e 

-  3.0% 

f 

-  U.c# 

NTU 

-l8-22j* 

Nst 

-l5-22j£ 

N3tNPr2/3 

±18~22?S 

NTU 

Uncertainty  at  20:1  odds 

6 

»22$ 

10 

«20# 

16 

-18% 

20 

Zl9fi 

1U 


DISCUSSION  OF  RESULTS 

To  evaluate  bh   >erformance  of  th*  test  section, 
60  x  60  x  0.0070  and  10  x  10  x  0,025  18/8  stainless  steel 
screen  matrici  .  ihic      simi   r  to  those  used  by  other 
inv  stigators  were  tested  (see  Table  III   and  Fig.  9)»      Fig.  1 
compares  the  results  of  the  present  investigat:  >n  with  th     '. 
interpretation  of  the  data  obtained  by  Coppage  and  Tong  (9) 
also  on  10  x  10  x  0„025)  screen  matrices.  Fig.  2  compares  the 
results  of  the  present  investigation  of  60  x  60  x  0.0070 
screens  with  the  best  interpretation  of  the  data  obtained  by 
Coppag   id  Tong  on  60  x  60  v  0.0075  scr  ns,  Th  heat 
transfer  and  fluid  friction  characteristics  obtained  for  the 
screen  matrices  are  well  /ithin  the  range  of  uncertainty  of 
reproducibility  of  Coppage,  Ton,  and  the  author. 

In  addition  to  the  screen  matrices,  which  wer  tested 
primarily  to  evaluate  the  performance  of   the  t    section, 
three  different  geometry  matrices  were  investigated:  two  of 
the  plate-fin  type  (Fig.  10)  and  one  fin  type  (Fig,  11)  in 
which  the  elements  were  skewed  to  prevent  meshing  of  "he  fins. 
Table  VIII  lists  the  details  of  these  matrices, 

Results  for  the  Type  U30  stairless  steel  plate-fin 
matrix  are  shown  in  Fig.  3  and  listed  in  Tables  III  and  IV, 
From  the  figure  it  is  s  n  that  in  the  higher  lieynolds  number 
range,  the  heat  transf  r   h  vior  follows  the  characteristic 
straight  lin<  tfith  a.  slope  of  slightly  less  than  the  minus  one 
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expectr  i  for  laminar  flow  in  smool    mnels.  it  a  ieynolds 
number  of  about  100  and  a  test  NTU  of  about  li>,  the  data  begins 
to  depart  I    the  straight  line  behavior.  As  the  Reynolds  num- 
ber decreases,  this  departure  becomes  more  pronounced  and  finally 
reaches  a  point  tfhere  Lt  appears  thai  the  heat  transfer  coef- 
ficient becomes  constant  or  independent  of  the  Ueynolds  number • 
This  behavior  is  due  to  the  longitudinal,  heat  cone  action  within 
the  matrix,  and,  as  mentioned  earlier,  was  expected  to  start 
showing  up  when  the  NTU  theoretically  had  a  value  of  11.7?(7). 

Since  the  conduction  parameter"""  for  a  given  geometry 
and  material  varies  inversely  as  the  Reynolds  number,  the 
smaller  the  Reynolds  number  the  greater  the  effect,   The  con- 
duction parameter  also  varies  inversely  as  the  flow  length, so 
that  a  first  observation  would  indicate  that  decreasing  the 
length  does  not  improve  the  situation;  but,  since  the  NTU  is 
directly  proportional  to  the  heat  transfer  area,  decreasing  the 
matrix  length  sufficiently  will  decrease  the  NTU  to  allow  test- 
ing a  lower  Reynolds  number  range.  In  this  investigation  the 
matrix  length  was  not  decreased  because  it  was  only  desired  to 
establish  the  straight  line  laminar  flow  behavior  and  sufficient 
lata  was  obtained  to  do  this  without  altering  the  matrix  length. 
Having  established  the  heat  transfer  behavior  in  the  higher  laminar 
Reynolds  number  range  thus  allows  extrapolation  of  the  data  into  ' 
low     ynolds  number  range,  it  should  be  noted,  however,  that  when 


conduction  parameter  is  defi  i*  i  as  .^2s_ 

wfCfL 
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such  information  as  this  is  utilized  in  heat  exchanger  design, 
conduction  plays  an  important  part  in  the  heat  exchanger's 
performance,  and  there  has  been  no  attempt  in  this  investigation 
to  obtain  correlations  which  ,/ould  predict  the  effects  of  con~ 
duction. 

The  friction  factor  curve  for  the  Type  U30  stainless 
steel  plate  fin  matrix  shows  essentially  a  straight  line 
behavior  with  a  slope  approximately  that  of  the  heat  transfer 
curve  as  would  be  expected  for  laminar  flow  conditions. 

Results  for  the  aluminum  plate-fin  matrix  are  shown 
in  Fig.  h   and  Table  V.  For  this  matrix  the  only  geometry  change 
was  that  of  the  hydraulic  diameter,  which  was  accomplished  by 
reducing  the  fin  height  and  pitch.  This  change,  along  with  the 
use  of  thinner  material,  provided  for  a  greatly  increased 
ratio  of  heat  transfer  area  to  unit  volume  and  of  free  flow  area 
to  frontal  area  compared  with  the  Type  U30  stainless  steel  plate- 
fin  matrixo  The  experimental  results  show  the  same  type  behavior 
as  previously  discussed  with  the  conduction  effect  beginning  at 
a  Reynolds  number  of  about  200  where  the  test  NTU  was  about  IJg 

In  Fig.  £  and  Table  VI  are  shown  the  results  for  the 
more  radical  of  the  geometry  changes.  For  this  matrix  the  plates 
were  removed  leaving  only  the  finsj  and  to  prevent  meshing  of 
the  fins,  the  elements  were  made  so  that  the  total  angle  between 
two  elements1  flow  channels  was  2!?°.  This  angle  was  chosen 
simply  because  it  was  an  easy  on=  to  fabricate  and  sufficiently 
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removed  from  a  zero  angle  to  produce  a  large  effect.   It  was 
believed  that  with  such  an  obviously  tortuous  f] ow  channel  the 
shape  of  the  result  curves  would  more  nearly  approach  those  of 
the  screen  matrices.   That  this  was  not  the  case  is  readily 
evident  by  comparing  Fig.  5   with  either  Fig.  1  or  2.   Instead, 
the  behavior  of  this  matrix  .•oras  quite  similar  to  that  of  the 
plate-fin  type  and  will  be  discussed  in  more  detail  later.  The 
conduction  effect  is  seen  to  begin  at  a  Reynolds  number  of  about 
3C0  for  which  the  test  NTU  was  about  12 „ 

For  comparison  purposes  the  best  interpretation  of  the 
heat  transfer  characteristics  of  the  three  finned  matrices  have 
been  replotted  and  extended  in  Fig.  6,  and  the  fluid  friction 
characteristics  reproduced  in  Fig,  7<? 

Comparing  the  stainless  steel  plate-fin  and  aluminum 
plate-fin  matrices  shows  the  effect  of  reducing  the  hydraulic 
diameter  by  reducing  size  and  pitch.  The  smaller  finned 
aluminum  matrix  has  approximately  the  same  fluid  friction 
characteristics  as  the  stainless  ^teel  matrix,  while  the  heat 
transfer  characteristic-  are  raised  hl%»     This  would  indicate 
that  the  fins  in  plate-fin  heat  exchanger?  could  be  reduced  in 
size  to  the  limit  allowed  by  considerations  other  than  heat 
transfer. 

The  comparison  of  the  aluminum  plate-fin  and  skewed  fin 
matrices  of  the  same  fin  size  shows  the  effect  of  removing  the 
plate  and  skewing  by  2£°,  These  changes  increase  the  fluid 
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friction  characteristics  by  100$  and  the  heat  transf 
characteristics  by  73o!?a»  Although  the  heat  transfer  rate  is 
increased,  this  may  be  an  expensive  means,  especially  in  gas 
turbine  regenerators  where   ]  '   increase  in  pressure  drop  due 
to  flow  friction  costs  more  than  a  \%   increase  in  specific  fuel 
consumption.  Possibly  there  is  a  skew  angle  at  which  the  in- 
crease in  heat  transfer  characteristics  is  sufficiently  greater 
than  the  increase  in  friction  characteristics  to  warrant  its 
use. 

Comparing  the  stainless  steel  plate-fin  and  skewed 
aluminum  fin  matrices  shows  an  increase  of  15>°^  in  the  heat 
transfer  characteristics  and  an  increase  of  as  much  as  120??  in 
the  fluid  friction  characteristics.  Ln   this  case,  the  heat 
transfer  area  per  unit  volume  is  nearly  the  same  for  these 
surfaces,  but  the  hydraulic  diameter  for  the  skewed  fins  is 
practically  5>0,b  greater,  indicating  that  there  would  be  less 
chance  of  fouling  with  the  skewed  arrangement. 
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CONCLUSIONS 

1.  A  successful  test  section  has  been  designed  and  fabricated 
to  obtain  connective  heat  transfer  and  friction  coefficients 
for  small  modules. 

2.  Connective  heat  transfer  end  friction  coefficients  have  been 
obtained  for  some  screen,  triangular  fin  and  plate-fin 
matrices  of  various  geometries. 

3 e  The  experimental  results  of  heat  transfer  behavior  and 
friction  factors  for  the  plate-fin  and  fin  matrices, 
summarized  in  Figs.  '6  and  7>  indicate  that  there  may  be 
better  geometries  for  compact  heat  exchanger  applications 
than  plate-fin  types  for  comparable  ratios  of  heat  transfer 
area  to  unit  volume. 
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Pin  pitch  =  19,8  per  lnoh 

Pin  height  =  ,0309  In, 

Fin  metal  thlokneea  ■  .002  in. 

Plate  metal  thlokness  ■  ,004  in, 

Matrix  In  direction  of  flow  =  3  in. 

Plow  pasaape  hydraulic  diameter  =  ,00177  ft,      „  « 

H«at  transfer  area/Unit  volume  of  matrix  -   1788  ft  ft 

Free  flow  area/prontal  area  =  ,7912 


^  -  present  data 
best  Interpretation 
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K  N  2/3  and  t   versus  NR.  for  430  stainless  steel  plate-fin 
Figure  5 
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i       NDIX  1 
SAMPLE  CALCULi     I  RUN  #1,     J  LATE- 

itecorded  Data 

tfi  ■  7Q«0°F.    (  u  rature) 

APorifice  *  I*. 62  in.  H?0 

^matrix  =  10.2  in.   H20 

^a  s  JO. 21  in.  Hga  (atmospheric  pressure) 

d  ■  2.310  in.  (orifice  diameter) 

D  s   3.03  in.  (pip  '        ':    c) 

t^-tfi  *  2.73  mv  (initial  temperature  difference,  5  thermocouples) 

chart  speed  a  2.0  sec/in 

chart  scale  -   0.558  mv/ir 

max.  slope  -  -6.772  in/in 

w*s  *  G«33k  lhm  (matrix  mass)  (Table  VIII) 

cs  -  0.216  3TU/lbm  °b\    (matrix  specific  heat)  (Table  YLLI) 

Dh  -  0.00.3^  ft.  (matrix  hydraulic  diameter)  (Table  VIII) 

Ac  »  0.051 91*  ft2  (matrix  flow  area)  (Table  VIII) 

a  ■  1*3.1*0  ft2  (matrix:  h  a1  transfer  area)  (Table  VIII ) 
Flow  Itate  Determination 

Using  the  method  in  (3)>  the  fluid  flow  rate  will  be  found 

wf  -  359  K  d2  Fa  Y  Y/°A?orificel  lbm/hr   I  (3),  p.  57,  equation  3 J 

K  =  flow  coeff     :,  velocity  of  approach  factor  includ 
(dimer>  i  nless) 

Fa  -  thermal  expansion  factor  (dimensionless) 

X  s  expansion  factor  (dimensionless) 

Ui 


/£>  *  density  of   I  I  !  i  it  in]  <       »rifice   (lbm/ft^) 

3  Rey  number  based  on  imeter 

/9  -   i/D   *  0.7^ 
Assuming  R-p  -  9i|,000 

K  -  0.7)473,         |(3),   p.   26,   Table  Sj 

d2  -  5.336   in2 

I  -  0.99*9       [(3),  p.  71,   Fig.  UObJ 


JJl.  -  ^o.ajOf9l)(/44)    _  676-^  /^ 


/°-°^w^   ^--7^)c^V 


Fa  -  1.000       |p),  ft.  67,  Fig.   33J 


wf  -  (3S9Ko.7U73)(S.336)(l.OOO)(0.99£>)       /  (0.0756K4.62) 
'    «  839  lbm/hr 

RD  -  0.O0U2UU  Wf/D/»       |("3)   p.  53,  equation  8J 

/4«  1.225  x  10"^  lbm/ft-sec,   at  70  °F,  [(10),   p.   69 

Table  2-8J  L 

aD  -  (0.00U2I1U) (839 )/(3. 08) (1.225  x  10"^)  =  9U,326 
The  value  of  R  *  9U,326  does  not  change  the  initial  flow     'icient, 

iyil'u  i.^U.r.-u i  na t  i  ..:. 

ro  obtain  NTU,  it  is  necessary  to  determine  the  experi- 
mental ma      .    of  the     •  Li   3  slope  of  the  generalized 
heatin,     ire: 


A(T/HTO) 


tA  AX 


1*2 


From  the  d 


Vv£  cs 


9  - 


and 


NJU  ~- 


h  A 

UTf.  Cf.     > 


'  . 


1*™  ■  ^7  (e  -  ^) 


and 


dil^NTu)--   fy±L\Jle 


ill  SO 


+  i 


and 


•«(IS=lf)-  <*nfe>*  <*»-'«>. 


ihus  the  expression  for  the  maximum  slope  becomes 


c£ 


.fcfl-g 


^tXntu) 


^H-tn) 


^  a 


JMAK 


MA* 


re  the  quantity 


MS* 


U3 


Lope  o  I  ture  d  Lf~ 

.   rsus  titne  curve.     As  this  t         r   I  tr  is 

13  a  f  Line,   the  i  1  >p 

directly  from  th      i  se.  Since  the  slope  is 

measured  directly  fro      hart,    th     pu1  tions  it  is 

necessary  to  introduce  tl    i  rt  scale— factor,  , 

cf  "   0.2)40h  BTU/lbm  °F,  fluid  specific  heat  at  70  °F, 
[(10),  p.  );1,  Table  2-3] 

-   /.  £6V  sec 


JUt^NTu) 


./-fsw)^4-77^^];0-8^ 


NTu  -  9.165  (Table  IK  -  Appendix  II) 

Avera^.-:  fluid  temperature  determination 

tfl  *  70  <V.  -  1.07  mv 

t:]_-tfi  a  2.73  ^v  (£  thermocouples) 

H-t.fi  ■  iillH  *  2ft?3  -  0.5U6  mv'(l  thermocouple) 

ti  "(ti-tfj+tfi  -  0.5U6  +  1,07  -  1.616  mv  -  88e87  1 

*ivg.  -  V'-fl  -  70  ■►  88.87  =  97.kU°F. 
2     -2 

Heat  transfer  Parameter  Determination 

Nst  -  NTU  fc  -  (9«l6jQ (0.059k)  =  o,0126 
a      TTToTio   " 

2/3  -  0.??5 ,  at  tavg,    ft  10),  p.  fl,  Table  2-lo] 
NstNPr2/3  -  (0.0126)1  ,   )  »  0.00958 


I* 


..  ,;.)'■> Ids  Nuitji' 

> 

N   -  Dl 
7T* 


.  •- 


G  -  "I 


(339) (3600 )/(0.059li  -T  3.92  lbm/sec-ft2 


/T  -  1.225  x  10-5  lbm/ft-sec,  at  t_v„    [(10),  p.  69, 
Table  2-8] 

■v.-   .  135K3.92)  /  (1.225  x  10-?)  =  , 
Friction  Factor  Determination 


P&M% 


rhe  friction  data  were  taken  at  anbient  a:  \  - 

/**-JL-  !■»  ■ 

v.f     R  Tfl 

*ave  ■  pl  -  ^  matrix 
2" 

pavg  »   (30.:  : .)(13.6)   -  10.2/2  =  1+05.9  in.  H20 

Pi  -  (U05.9)(5.20U)/(53.3)(530)  -   907lt?  lbm/ft-3 

f  »  (lO.2)(5.2OU)(O.O7U7)(2)(32.2)(O.059t)"/  (3.92)2( 
=    0.0225 


U5 


.  £  ..IX  II 
;   ;.  I  ..  [Ml      L    .         "  .  .     E     :ng  cu.,. 

Lock,    (l)   derived  the  e;-  Lon  for  th:    slope  of  the 

generaliz  i  irvc  eqv       o  page  ?)>''- 

X    *    L    t.;-:j.t    _.    -    TvTU, 


" 


_____ 


~ 


^(       ,■       ■ 


ly/NT  ■    T 


. 


i,  IDC  JL60lj  digital  computer  was    ised  tc        L       I  ■.■• 


.  \ 


•  : ,'  :  ■■.    I  by 


v 
of       '  i  _        _r 

_  !    / 

varying  f/lSTU   for  particular  values  of  NTU.  For        r  NTU  -  2, 
and  using  an  increm  nit  of        :  .00001,     i  xi  ium  sloj: 


...      he  initial     ■■  ■■  r       /NTU  ei  ■    Ly 


. 

.-  NTU  >  2,  incr    its  of  '      ■  .00]  ■  used.      /NTU 

Ltil  a  val*   of  tl    lope  sms  .  than  the  pn  v_  >us 

value  was  obtained,  and  this  maximum  slope  was  recorded.  The  com™ 

puter  then  moved  to  the  next  value  of  MTU, 

1______1- 

The  values  of  1    ~  for  values  of  NTU 

and   f  •  /NTU  are  cor.ta:*: 


Table U 
Values   of    the   Max.    Slope    of    the    Generalized   Heating  Curve  ,NTU  and  yl/WU 


Qfth-*^ 

w 

NTU 

TNTU 

tea 

Ml 

NTU 

Ajtu 

Ml 

'NTU 

CjjzQ 

«h 

NT1J 

/ntu 

It,,  - t. ) 

1   NTU 

[dC^WTvJ 

[dCVuTu). 

Ld(^*/TU,],J 

[<ny*Tu\ 

.00009V 

.0  1 

.      1          1 

.19318,' 

1  .  u  "> 

.0000 1 

.571957 

2.91 

,*M020( 

.60  3012 

1.37 

.62800 

.000)9/ 

.0? 

.1)000  1 

.190031 

I.U6 

■  1 

.572166 

2.V2 

.10101 

.651601 

1  .  )8 

.62900 

.0008/1 

.0-1 

.   0001 

.196811 

I.U7 

.00001 

.572977 

2.7  3 

.M070O 

.601166 

1.19 

.6  5000 

.001537 

•UOOOI 

.198616 

I.UH 

.uooo.  1 

.5/3189 

2.  71 

.10900 

.601/28 

1.10 

.65100 

.0023/6 

.05 

.5000  1 

.000318 

I.U9 

.uOOO 1 

.5/1005 

2.90 

.-*  1200 

.600290 

1.11 

.65200 

.005390 

.06 

.UOOOI 

.002011 

I.SO 

.00001 

.5/1518 

2.96 

.1*  UOO 

.605851 

1.12 

.6  5500 

.001067 

.07 

.  i^OOO  1 

. 00  1690 

l.si 

.UOOOI 

.5/0050 

. 

.-*  1700 

.6061 12 

1.13 

.63100 

.000906 

.0  \ 

.uCOOl 

.000510 

1.52 

.0000  1 

.070001 

.4 1900 

.606971 

1.11 

.63500 

.00/105 

.09 

.'«000  1 

.506PB7 

1.5  3 

.UOOOI 

.5/60/3 

.    .    ■  1    :■ 

.60/031 

1.10 

.6  5600 

.00/01. 

.  1(1 

.',U(  (!   1 

.008121 

1.51 

.UOOOl 

.5/6091 

1 .  00 

.42100 

.658095 

.63/00 

.010610 

.  1 1 

.uOOO 1 

.009921 

1.55 

.00001 

.0/71 1 1 

.12600 

.658656 

1.17 

.65800 

.0  12//.' 

.  12 

.00001 

.01 1580 

1.56 

.00001 

.0//61 1 

1.0  2 

.42900 

.609216 

1.18 

.65900 

.0  I18ii 

.  1  1 

.     00  0  1 

.012809 

1.S7 

.078166 

5.03 

.15100 

.609776 

1.19 

.65900 

.017034 

.uOOO 1 

.0 11 190 

1.58 

.OCOOI 

.578691 

3.01 

.4  3300 

.660336 

1.00 

.61000 

.0IV566 

.uooo 1 

.0  10015 

1.5  t 

.UOOOI 

.079220 

.4  3S0O 

.660896 

1.0  1 

.61100 

.0/ It  16 

•UOOO 1 

.0  1680) 

.0000 1 

.079719 

5.O6 

.661155 

1.02 

.61200 

.021.  id/ 

'.  1  7 

.!    OOOI 

.018127 

1.6  1 

.0000  1 

.580279 

1 . 0  l 

.1*4000 

.662011 

1.0  1 

.61300 

.02/065 

.0000 1 

.019361 

1.62 

.00001 

.58081 1 

.44200 

.662073 

1.51 

.61100 

.02/655 

•UOOOl 

.020061 

1  .'.  S 

.UOOOI 

.58  151) 

5.09 

.44100 

.661132 

1.55 

.61500 

.012/1V 

.UOOOI 

.02  1727 

1  .Oil 

.oOOOl 

.0818/6 

3.  10 

.1470o 

.66  1691 

1.06 

.61600 

.01WX 

'.2  1 

.1  0001 

.022801 

1  .6S 

.uooo 1 

.08/1 1 1 

3.  1  1 

.4490C 

.661219 

U.07 

.61600 

.056012 

.22 

.0000 1 

.S25910 

1  .  66 

.00001 

.582917 

1.  12 

.45100 

.66180/ 

1.08 

.61/00 

.012051 

.23 

.UOOOl 

.020001 

1.67 

. 

.58  516) 

5.  1  1 

.4  5300 

.665365 

1.09 

.61800 

.01511, 

. .  4 

.UOOOI 

.026020 

.JOOO 1 

. 58*02 1 

3.  11 

.45500 

.660923 

1.60 

.61900 

.0106  76 

.25 

.oOOOl 

.02  7000 

1  .6V 

.000 0  i 

.581060 

3.  15 

.4570C 

.666180 

1.6  1 

.60000 

.002121 

.26 

.uooo 1 

.027901 

1.  70 

. ouoo  1 

.580099 

3.  16 

.45  /(' 

.667038 

1.62 

.65100 

.000600 

.27 

.0000 1 

.028868 

1.71 

1 1 

.585610 

5.17 

.46  100 

.66/595 

1.61 

.651011 

.00920  1 

.7' 

.oOOOl 

.029718 

\.l? 

1 

.586  18  1 

1.18 

.46401 

.668101 

1.61 

.60200 

.06292/ 

.2  i 

.  jOOO 1 

.050591 

\  .1  '. 

.0000 1 

.086/2) 

i.W 

.   .  ■  - 

.668/08 

1.60 

.60300 

.0666/1 

.  so 

.UOOOI 

.01 1101 

1  .  fU 

. 

.08/26/ 

.46^00 

.669261 

1.66 

.65100 

.0/0161 

.  3  1 

.0000 1 

.012182 

1.  r5 

.00001 

.08/8 1 1 

1.2  1 

.47000 

.66/8/  1 

1.67 

.6550C 

.  o/i 16/ 

1  1 

.032920 

I  .  ft, 

.088)00 

1.22 

■ 

.6/0376 

1.68 

.60600 

.0/829 I 

'.  .  s 

.'000  1 

.015650 

1.77 

1.21 

.6/09  52 

1.69 

.60600 

.082760 

.uooo I 

.051512 

1  .  7* 

.089118 

3.21 

1 17500 

.6/118/ 
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.60/00 

.066121 

'.  35 

.00001 

.0  11906 
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...  7  700 

.6/2015 
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